Acoustic radiation- and streaming-induced microparticle velocities determined by 

micro-PIV in an ultrasound symmetry plane 



Rune Barnkob( a )*, Per Augustsson( b ' , Thomas Laurell^ b,c ^, and Hcnrik Bruus ( - d ) 
^Department of Micro- and Nanotechnology, Technical University of Denmark 

DTI] Nanotech Building 345 East, DK-2800 Kongens Lyngby, Denmark 
^Department of Measurement Technology and Industrial Electrical Engineering, 
Division of Nanobiotechnology, Lund University, Box 118, S-221 00 Lund, Sweden 
^Department of Biomedical Engineering, Dongguk University, Seoul, South Korea 
^Department of Physics, Technical University of Denmark 
DTU Physics Building 309, DK-2800 Kongens Lyngby, Denmark 
(Dated: 31 August 2012) 

We present micro-PIV measurements of suspended microparticles of diameters from 0.6 um to 
10 jim undergoing acoustophoresis in an ultrasound symmetry plane in a microchannel. The motion 
of the smallest particles are dominated by the Stokes drag from the induced acoustic streaming 
flow, while the motion of the largest particles are dominated by the acoustic radiation force. For 
all particle sizes we predict theoretically how much of the particle velocity is due to radiation 
and streaming, respectively. These predictions include corrections for particle-wall interactions 
and ultrasonic thermoviscous effects, and they match our measurements within the experimental 
uncertainty. Finally, we predict theoretically and confirm experimentally that the ratio between 
the acoustic radiation- and streaming-induced particle velocities is proportional to the square of 
the particle size, the actuation frequency and the acoustic contrast factor, while it is inversely 
proportional to the kinematic viscosity. 
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I. INTRODUCTION 

Acoustofluidics and ultrasound handling of particle 
suspensions, recently reviewed in Review of Modern 
Physics [l| and Lab on a Chip 0, is a field in rapid 
growth for its use in biological applications, such as sep- 
aration and manipulation of cells and bioparticles. Mi- 
croChannel acoustophoresis has largely been limited to 
manipulation of micrometer-sized particles, such as yeast 
3f, blood cells [||, cancer cells [EH3, natural killer cells 
8jj, and affinity ligand complexed microbeads 0], for 
which the acoustic radiation force dominates. Precise 
acoustic control of sub-micrometer particles, e.g. small 
bacteria, vira, and large biomolecules remains a chal- 
lenge, due to induction of acoustic streaming of the sus- 
pending fluid. Nevertheless, acoustic streaming has been 
used to enhance the convective transport of substrate 
in a microenzyme reactor for improved efficiency [Io| . 
while acoustic manipulation of sub-micrometer particles 
has been achieved in a few specific cases including en- 
hanced biosensor readout of bacteria [ll| and bacterial 
spores [lH , and trapping of E-coli bacteria [HI . 

When a standing ultrasound wave is imposed in a mi- 
crochannel containing an aqueous suspension of particles, 
two forces of acoustic origin act on the particles: the 
Stokes drag force from the induced acoustic streaming 
and the acoustic radiation force from sound wave scatter- 
ing on the particles. To date, the experimental work on 
acoustophoresis has primarily dealt with cases where the 
acoustic radiation force dominates the motion, typically 
for particles of diameters larger than 2 um. Quantitative 



experiments of 5-pm-diamcter polymer particles in wa- 
ter [I3 - [l6l | have shown good ag reement with the classical 
theoretical predictions [13, 13 °f the acoustic radiation 
force acting on a microparticle of radius a much smaller 
than the acoustic wavelength A, and where the viscos- 
ity of the suspending fluid is neglected. However, as the 
particle diameter is decreased below 2 um, a few times 
the acoustic boundary-layer thickness, the particle mo- 
tion is typically strongly influenced by the Stokes drag 
force from the induced acoustic strea ming flow, which 
has been reported by several groups [l9l - [2l1 | , and the ra- 
diation force is modified due to the acoustic boundary 
layer [H||. 

As pointed out in a recent review [jjjj . the acoustic 
streaming is difficult to fully characterize due to its many 
driving mechanisms and forms. In acoustofluidic sys- 
tems, the streaming is primarily boundary-driven arising 
at rigid walls from the large viscous stresses inside the 
sub-micrometer-thin acoustic boundary layer of width 
5. The boundary-driven acoustic streaming was theo- 
retically treated by Rayleigh [23| for an isothermal fluid 
in an infinite parallel-plate channel with a standing wave 
parallel to the plates of wavelength A much larger than 
the plate distance h, and where h is large compared 
to 8, i.e. X h 5. However, in many applications 
of acoustofluidic systems the channels provide enhanced 
confinement, the acoustic wavelength is comparable to 
the channel height, and the liquid cannot be treated as 
being isothermal. Rayleigh's prediction is often cited, but 
to our knowledge the literature contains no quantitative 
validation of its accuracy when applied to acoustofludic 



2 



systems. This lack of quantitative tests is most likely 
due to the fact that boundary-driven acoustic streaming 
is very sensitive to geometry and boundary conditions, 
making it difficult to achieve sufficient experimental con- 
trol. However, quantitative comparisons between theory 
and experiment of acoustic streaming are crucial for the 
advancement of the acoustofluidics research field. Un- 
derstanding and controlling the ratio of radiation- and 
streaming-induced acoustophoretic velocities may be the 
key for future realization of ultrasound manipulation of 
sub-micrometer particles. 

In 2011, we presented a temperature-controlled 
micro-PIV setup for accurate measurements of the 
acoustophoretic microparticle motion in a plane (l5j |. 
Here, we use the same system and the ability to estab- 
lish a well-controlled transverse resonance for quantita- 
tive studies of how much the radiation- and streaming- 
induced velocities, respectively, contribute to the total 
acoustophoretic velocity. More specifically, as illustrated 
in Fig. [TJ we study the microparticle motion in the ul- 
trasound symmetry plane (magenta) of a straight rect- 
angular microchannel of width w = 377 pm and height 
h = 157 pm. We determine the velocities for parti- 
cles of diameter 2a ranging from 0.6 pm to 10 pm, and 
based on this we examine the validity of Raylcigh's the- 
oretical streaming prediction. We also derive theoreti- 
cally and validate experimentally an expression for the 
microparticle velocity as function of particle size, ultra- 
sound frequency, and mechanical properties of the sus- 
pending medium. 
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FIG. 1. (a) Sketch of the silicon/ glass microchip used in 
our experiments, see also Ref. It contains a straight 

rectangular water-filled microchannel (light cyan) of length 
L = 35 mm, width w = 377 pm, and height h = 157 pm. 
(b) The horizontal ultrasound symmetry plane (magenta) of 
length I = 892 pm and width w in the xy plane at the center 
of the channel, (c) The vertical channel cross section (cyan). 



streaming flow. Both effects arc time-averaged second- 
order effects arising from products of the first-order fields. 
The drag force from the acoustic streaming flow domi- 
nates the motion of small particles, while the motion of 
larger particles are dominated by the acoustic radiation 
force. This is clearly illustrated in recent numerical sim- 
ulations by Muller et al. 2J|, which are reproduced in 



Fig. [2] (b) the streaming flow advects small particles in 
a vortex pattern, and (c) radiation force pushes larger 
particles to the pressure nodal plane at y = 0. 



II. THEORY OF SINGLE-PARTICLE 
ACOUSTOPHORESIS 



In this work we study a silicon-glass chip containing 
a rectangular microchannel sketched in Fig. [T] and de- 
scribed further in Section IIIII The microchannel con- 
tains a particle suspension, and the chip is ultrasoni- 
cally actuated by attaching a piczo transducer to the 
chip and driving it with the voltage U pp at the angu- 
lar frequency w — 2nf, where / is a frequency in the low 
MHz range. By proper tuning of the applied frequency, 
the actuation induces a resonant time-harmonic ultra- 
sonic pressure field pi(r) exp (— iuit) and velocity field 
Vi(r) exp (— kjt), here expressed in the complex time- 
harmonic notation. Throughout the paper, we only study 
the case of a ID transverse pressure resonance of ampli- 
tude p a and wavenumber k = 2ir/X, 



Pi( r ) = V a 



cos 



k[y- 



(1) 



The case of A/2 = w or k = tt/iv is shown in Fig. \^&). 

The particle suspensions are dilute enough that the 
particle-particle interactions are negligible, and thus only 
single-particle effects are relevant. These comprise the 
acoustic radiation force due to particle-wave scattering 
and the viscous Stokes drag force from the acoustic 



The acoustic radiation force 



We consider a spherical particle of radius a, density 
p p , and compressibility K p suspended in a liquid of den- 
sity p , compressibility k , viscosity 77, and momentum 
diffusivity v = T]/p . Recently, Bruus and Settnes [22[ 
gave an analytical expression for the viscosity-dependent 
time-averaged radiation force _F rad in the experimentally 
relevant limit of the wavelength A being much larger than 
both the particle radius a and the thickness 8 = \j2v /uj 
of the acoustic boundary layer, without any restrictions 
on the ratio 8 = 8 /a. For the case of a ID transverse pres- 
sure resonance, Eq. (TT]), the viscosity-dependent acoustic 
radiation force on a particle located at (y, z) reduces to 
the z-independent expression 



(y, z) = Wfc£ ac p, 8) sin (2ky) , (2) 



where E a 



jn pi is the time- averaged acoustic energy 
density and where the acoustic contrast factor $ is given 
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FIG. 2. Numerical simulation of microparticle acoustophoresis in the vertical microchannel cross section of Fig. djc) adapted 
from Ref. [24| . The magenta line represents the ultrasound symmetry plane of Fig. \V[b) . (a) Color plot of the transverse standing 
ultrasound pressure wave pi of Eq. |TJ) ranging from — p a (dark blue) to p a (dark red), (b) Trajectories (blue lines) of small 
0.5-um-diameter particles (dots) dominated by the Stokes drag force from the boundary-induced streaming, (c) Trajectories 
(colored lines) of large 5.0-um-diameter particles (dots) dominated by the acoustic radiation force. 



in terms of the material parameters as 
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We note that for all the microparticle suspensions studied 
in this work including the viscous 0.75:0.25 water: glycerol 
mixture, the viscous corrections to $ are negligible as we 
find | 5)/$(k,£,0) 1 1 < 0.4 %. 

If F Tad is the only force acting on a suspended parti- 
cle, the terminal speed of the particle is ideally given by 
the Stokes drag as u rad = F lad / (6-KTja) . Using Eq. © 
for the transverse resonance, M rad only has a horizontal 
component w5, ad , and this can be written in the form 



rad 



*0 2 



sin(2fcy), 



(4) 



where the characteristic velocity amplitude u and par- 
ticle radius a are given by 
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Here, Zq is the characteristic acoustic impedance, and 
the numerical values are calculated for polystyrene par- 
ticles suspended in water using parameter values listed 



in Section HU with 
Barnkob et al. \U 



2 MHz and E a 



10 J m 6 as in 



B. Boundary-driven acoustic streaming 

In 1884 Lord Rayleigh [23| published his now classical 
analysis of the boundary-driven acoustic streaming veloc- 
ity field U in an infinite parallel-plate channel induced 



by a first-order bulk velocity field having only a horizon- 
tal y-component given by v\ = v a sin [k(y-w/2)]. This 
corresponds to the first-order pressure of Eq. {T|) illus- 
trated in Fig. [UJa). For an isothermal fluid in the case of 
A ^> h ^> S, Rayleigh found the components U y and U z 
of U outside the acoustic boundary to be 
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A plot of U driven by the ID transverse standing half- 
wave resonance is shown in Fig. [3] We expect this ana- 
lytical expression to deviate from our measurements be- 
cause the actual channel does have side walls, it is not 
isothermal, and instead of A > /i we have A = 4.8ft for 
A = 2w and A = 2Ah for A = w. 

At the ultrasound symmetry plane z = 0, U only has a 
horizontal component, which we denote u^ tr . In analogy 
with Eq. this can be written as 



uf = u s° p sm(2fcy), 



*p = Ye * °' 188 ' 



(7) 



where the sub- and superscript in the streaming coeffi- 
cient s° refer respectively to the parallel-plate geometry 
and the isothermal liquid in Rayleigh's analysis. 

To estimate the effect on u^ tr of the side walls and the 
large height h ps A in the rectangular channel of Fig. [IJc) , 
we use the numerical scheme developed by Muller et al. 
[HI for calculating the acoustic streaming based directly 
on the hydrodynamic equations and resolving the acous- 
tic boundary layers, but without taking thermoviscous 
effects fully into account. The result shown in Fig. 0] re- 
veals that u S y T is suppressed by a factor of 0.82 in the 
rectangular geometry relative to the parallel-plate geom- 
etry and that it approaches zero faster near the side walls 
at y = ±ii)/2. The approximate result is 



< tr « u o s r sin(2fcy), 



0.154, 



(8) 
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FIG. 3. Vector plot (black arrows) of the acoustic streaming 
U and color plot of its magnitude U from zero (dark blue) to 
maximum (dark red) given by Rayleigh's analytical expression 
Eq. ((6| and valid for a shallow isothermal liquid slab (A ^> h) 
between two parallel plates (gray) of distance h driven by a 
standing ultrasound pressure wave p x — p a cos \k(y — w/2)] 
for A/2 = w or k = n/w with w = 377 um. 



where the sub- and superscript in the streaming coeffi- 
cient s° refer respectively to the rectangular geometry 
and the isothermal liquid. 

We estimate the thermoviscous effect on u y tr , in partic- 
ular the temperature dependence of viscosity, using the 
analytical result by Rednikov and Sadhal for the parallel- 
plate geometry [251 ] . They found a streaming factor sp" 
enhanced relative to s^, 
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where a is the thermal expansion coefficient, D th the 
thermal diffusivity, and 7 the specific heat ratio, and 
where the value is calculated for water at T = 25 °C. 

Combining the reduction factor 0.82 from the rectan- 
gular geometry with the enhancement factor 1.26 from 
thermoviscous effects leads to sT 



1.03 sp 1 or 



sj sin(2fcy), 



0.194, 



(10) 



where the sub- and superscript in the streaming coeffi- 
cient s^T refer respectively to the rectangular geometry 
and a thermoviscous liquid, see Fig. |U 



C. Acoustophoretic particle velocity 

A single particle undergoing acoustophoresis is directly 
acted upon by the acoustic radiation force F Iad , while the 
acoustic streaming of velocity U contributes with a force 
on the particle through the viscous Stokes drag from the 
suspending liquid. Inertial effect can be neglected as the 
characteristic time scale p p a 2 jr\ of acceleration (< 1 ps) is 
minute in comparison with the time scale of the motion of 
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FIG. 4. Numerical simulation of the horizontal component 
Uy of the streaming velocity in the ultrasound symmetry 



plane at 



0. In this plane the vertical component is zero, 



u\ 1 = 0. Three cases are shown: the isothermal parallel-plate 
channel Eq. (0, the isothermal rectangular channel Eq. l[8]l. 
and the thermoviscous rectangular channel Eq. (|10|l . 



particles (> 1 ms). The equation of motion for a spherical 
particle of velocity u p then becomes 



prad 



U 



U. 



(11) 



As we have seen above, there are no vertical velocity 
components in the ultrasound symmetry plane at z = 0, 
and combining Eqs. (0} and ([TO")) we obtain the horizontal 
particle velocity component u p of amplitude w a , 



,rad 
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where we have dropped the sub- and superscripts of the 
streaming coefficient s. The ratio of the radiation- and 
streaming-induced velocity amplitudes becomes 
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which scales linearly with the angular frequency lo and 
the square a 2 of the particle radius, but inversely with 
the streaming coefficient s and the momentum diffusivity 
1//$ rescaled by the acoustic contrast factor. 

In Fig. [S] we show colored contour plots of the ratio 
u T y ad I 'u & y l ': in (a) for polystyrene particles in water at 
25 °C as function of the particle diameter 2a and the 
ultrasound frequency /, and in (b) as function of / and 
the rescaled momentum diffusivity f/$ for fixed parti- 
cle diameter 2a = 2 um. The green dots indicate the 
experiments described in Sections IIIII and IIVI 

We define the critical particle diameter 2a c for cross- 
over from radiation-dominated to streaming-dominated 
acoustophoresis as the particle diameter for which 



,rad 



(a c ) 



r (a c ). This results in 



2a c = 2ag y/s = 



2Av s 

ss 1.4 um, 



(14) 



where the numerical value is calculated for polystyrene 
particles in water ($ = 0.17) at / = 2 MHz using s = . 
For a = a c the ratio of the velocity amplitudes u™ d /u s }* is 
unity, and consequently the unity contour line in Fig.[5ja) 
represents 2a c as function of ultrasound frequency / 
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FIG. 5. Theoretical prediction of the ratio u™ d /u S y T Eq. (|13|l of radiation- and streaming-induced velocities in the ultrasound 
symmetry plane of the microchannel (magenta in Figs. Q] and [5J) at 25 °C. (a) Contour plot of uj, ad /«y tr as function of particle 
diameter 2a and ultrasound frequency / for a suspension of polystyrene particles PS in pure water, (b) Contour plot of «J / ad /«° tr 
as function of / and the rescaled momentum diffusivity z//$ for fixed particle diameter 2a = 2 p. The green dots indicate the 
cases studied experimentally in Sections IIIII and llVI Suspensions of polystyrene particles in three different liquids are indicated 
by dashed lines: pure water (light blue) as well as 0.75:0.25 and 0.67:0.33 mixtures of watenglycerol (white). 



TABLE I. The wall correction factor x to the single-particle 
drag for the particle sizes used in the experiment. 



2a 


paral 


paral 


pcrp 


0.6 pm 


1.004 


1.005 


1.004 


1.0 pm 


1.006 


1.008 


1.006 


1.9 pm 


1.012 


1.016 


1.011 


2.6 pm 


1.017 


1.022 


1.016 


4.9 pm 


1.032 


1.042 


1.030 


10.2 pm 


1.070 


1.092 


1.065 


D. 


Wall corrections to single-particle drag 



The sub-millimeter width and height of the rectangu- 
lar microchannel enhance the hydrodynamic drag on the 
microparticles. This problem was treated by Faxen for of 
a sphere moving parallel to a planar wall or in between 
a pair of parallel planar walls [26j and later extended by 
Brenner [27| to motion perpendicular to a sin gle planar 
wall, as summarized by Happel and Brenner [28| . The 
enhancement of the Stokes drag is characterized by a di- 
mensionless correction factor x( a ) modifying Eq. (1121) . 



1 



u sin(2fcy). 



(15) 



No general analytical form exists for x, so we list the 
result for three specific cases. For a particle moving par- 
allel to the surface in the symmetry plane z = in the 
gap of height h between two parallel planar walls, \ is 



paral 
Xz=0 



[l - 1.004(2a//i) + 0.418(2a//i) E 



1.070, 

(16) 



while for motion in the planes at z = ±h/4 it is 
[1 - 1.305(2a/fr) + 1.18(2a/fr) 3 ] _1 ft 



paral 
Xz = ±h/4 



s 1.092. 

(17) 

Here the numerical values refer to a particle with diam- 
eter 2a = 10 pm moving in a gap of height h = 157 pm. 
Similarly, for particle motion perpendicular to a single 
planar wall, the correction factor is 



X P CT P = -sinh(a)]r 



n(n + 1) 



3 ------ v , ^ (2n-l)(2n + 3) 

2 sinh [(2n + l)a] + (2n + 1) sinh (2a) 



4 sinh 2 [(n + |)a] - (2n + l) 2 sinh 2 (a) 



- 1 



1.065, 



(18) 



where a = cosh -1 (Ay /a) and Ay is the distance from the 
center of the particle to the wall. The numerical value 
refers to a 10- pm particle located at Ay = w/4. 

The values of the wall correction factor \ f° r ai l the 
particle sizes used in this work are summarized in TablcQ] 



III. EXPERIMENTAL PROCEDURE 

Experiments were carried out to test the validity of the 
theoretical predictions for the acoustophorctic particle 
velocity Eq. (|15[) in the horizontal ultrasound symmetry 
plane and for the ratio of the corresponding radiation and 
streaming-induced velocities, see Eq. (p~5|) and Fig. [5l 

We use the experimental technique and micro-PIV sys- 
tem as presented in Augustsson et al. The setup is 
automated and temperature controlled. This enables sta- 
ble and reproducible generation of acoustic resonances as 
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TABLE II. Material parameters at T = 25 °C. 



Polystyrene 



Density [29] 


Pp 


1050 kg m~ 3 


Speed of sound [30] (at 20 °C) 


c p 


2350 ms" 1 


Poisson's ratio [31J 


o-p 


0.35 


Compressibility^ 


K p 


249 TPa" 1 


Water 


Density [29] 


PO 


997 kgm -3 


Speed of sound [29] 


CO 


1497 ms" 1 


Viscosity [29] 


V 


0.890 mPas 


Viscous boundary layer, 1.940 MHz 


5 


0.38 urn 


Viscous boundary layer, 3.900 MHz 


S 


0.27 urn 


Compressibility^ 


K o 


448 TPa" 1 


Compressibility factor (polystyrene) 


h 


0.444 


Density factor (polystyrene) 


h 


0.034 


Contrast factor (polystyrene) 


$ 


0.17 


Rescaled momentum diffusivity (polyst.) 


vj<$> 


5.25 mm 2 s" 1 


0.75:0.25 mixture of water and 


glycerol 


Density [33] 


PO 


1063 kgm" 3 


Speed of sound [341 


CO 


1611 ms" 1 


Viscosity [33] 


V 


1.787 mPas 


Viscous boundary layer, 2.027 MHz 


6 


0.51 urn 


Compressibility^ 


K o 


363 TPa" 1 


Compressibility factor (polystyrene) 


h 


0.313 


Density factor (polystyrene) 


h 


-0.008 


Contrast factor (polystyrene) 




0.10 



Rescaled momentum diffusivity (polyst.) v/<& 16.8 mm s 1 
a Calculated as re D = — — —-, — ^tt from Ref. |32| . 

P l+CTp (p p C p ) 

b Calculated as k q = l/(p Cg) 

a function of temperature and frequency. It also enables 
repeated measurements that lead to good statistics in the 
micro-PIV analyses. The resulting acoustophoretic par- 
ticle velocities are thus of high precision and accuracy. 

Using the chip sketched in Fig. [TJ a total of 22 sets 
of repeated velocity measurement cycles were carried out 
on polystyrene particles of different diameters undergo- 
ing acoustophoresis in different suspending liquids and at 
different ultrasound frequencies. In the beginning of each 
measurement cycle, a particle suspension was infused in 
the channel while flushing out any previous suspensions. 
Subsequently, the flow was stopped, and a time lapse mi- 
croscope image sequence was recorded at the onset of the 
ultrasound. The cycle was then repeated. 

A. Microparticle suspensions 

Two types of microparticle suspensions were exam- 
ined; polystyrene particles suspended in Milli-Q water 
and polystyrene particles suspended in a 0.75:0.25 mix- 
ture of Milli-Q water and glycerol. To each of the two 



TABLE III. The nominal and the measured diameter of the 
polystyrene particles used in the experiment. 

Nominal diameter Measured diameter (2a) 

591 nm (0.59 ± 0.03) un£ 

992 nm (0.99 ± 0.05) un£ 

2.0 um (1.91 ± 0.07) un£ 

3.0 um (2.57 ± 0.07) un£ 

5 um (5.11 ±0.16) urn^ 

10 um (10.16 ± 0.20) un£ 

a Value from manufacturer and assumed 5 % standard deviation. 
b Measured by Coulter counter 

suspending liquids was added 0.01 % w/V Triton-X sur- 
factant. The material parameters of the suspensions are 
listed in Table HU Note that the rescaled momentum dif- 
fusivity v/$> of the glycerol suspension is 3 times larger 
than that for the Milli-Q water suspension. 

We analyzed 12 particle suspensions by adding parti- 
cles of 6 different diameters 2a from 0.6 um to 10 um 
to the two liquids. The particle diameters were mea- 
sured using a Coulter Counter (Multisizer 3, Beckman 
Coulter Inc., Fullerton, CA, USA) and fitting their dis- 
tributions to Gaussian distributions, see Supplemental 
Material. The resulting diameters are listed in Table Mil 

The concentration C of the particles were calculated 
based on the concentrations provided by the manufac- 
turer and varies in this work from 10 10 m -3 for the largest 
particles in the 0.75:0.25 mixture of water and glycerol 
to 10 15 m~ 3 for the smallest particles in the pure water 
solution. The concentrations correspond to mean inter 
particle distances C" 1 / 3 ranging from 4 particle diame- 
ters for the largest 10-um particle in water to 173 particle 
diameters for the smallest 0.6-um particle in the 0.75:0.25 
mixture of water and glycerol. Mikkelsen and Bruus [35| 
have reported that hydrodynamic effects become signif- 
icant for interparticlc distances below 2 particle diame- 
ters. Thus we can apply the single-particle theory pre- 
sented in Section |nj 



B. Measurement series 

We measured the acoustophoretic velocities of polysty- 
rene microparticles in the following four series of experi- 
ments, the second being a repeat of the first: 

MQO: Milli-Q water, / = 1.940 MHz, 

A = 2w, and 2a = 1.0, 1.9, 2.6, and 5.1 um. 
MQ1: Milli-Q water, / = 1.940 MHz, 

A = 2w, and 2a = 0.6, 1.0, 1.9, 2.6, 5.1, and 10.2 um. 
MQ2: Milli-Q water, / = 3.900 MHz 

A = w, and 2a = 0.6, 1.0, 1.9, 2.6, 5.1, and 10.2 um. 
GL2: 0.75:0.25 Milli-Q water:glycerol, / = 2.027 MHz 

A = 2w, and 2a = 0.6, 1.0, 1.9, 2.6, 5.1, and 10.2 um. 

Given the different particle diameters, we thus have 
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the above-mentioned 22 sets of acoustophoretic particle- 
velocity measurements, each consisting of 50 to 250 mea- 
surement cycles. All experiments were carried out at 
a fixed temperature of 25 °C and the applied piezo 
voltage U*. The camera frame rate was chosen such 
that the particles would move at least a particle diame- 
ter between two consecutive images. The measurement 
field of view was 1280 x 640 pixels corresponding to 
892 um x 446 jim. The imaging parameters were: optical 
wavelength 520 nm for which the microscope objective is 
most sensitive, numerical aperture 0.4, and magnification 
20. See acquisition details in the Supplemental Material. 



C. Micro-PIV analysis 

The micro-PIV analyses were carried out using the 
software EDPIV - Evaluation Software for Digital Par- 
ticle Image Velocimetry, including the image procedure, 
the averaging in correlation space, and the window shift- 



ing described in detail in Ref. [15[. For the MQ1, MQ2, 
and Gil series, the interrogation window size was 32 x 
32 pixels with a 50 % overlap resulting in a 79 x 39 square 
grid with 16 pixels between each grid point. For the MQ0 
series, the interrogation window size was 64 x 64 pixels 
with a 50 % overlap resulting in a 39 x 19 square grid 
with 32 pixels between each grid point. 

In micro-PIV all particles in the volume are illumi- 
nated, and the thickness of the measurement plane is 
therefore related to focal depth of the microscope objec- 
tive. This thickness, denoted the depth of correlation 
(DOC), is defined as twice the distance from the mea- 
surement plane to the nearest plane for which the par- 
ticles are sufficiently defocused such that it no longer 
contributes significantly to the cross-correlation analysis 
[36| . The first analytical expression for the DOC was 
derived by Olsen and Adrian (37l | and later improved by 
Rossi et al. (38|. Using the latter, we found that the DOC 
ranges from 14 um to 94 um sa h/2 for the smallest and 
the largest particles, respectively. Consequently, in the 
vertical direction all observed particles reside within the 
middle half of the channel. 



IV. RESULTS 

The core of our results is the 22 discrete acoustopho- 
retic particle-velocity fields obtained by micro-PIV anal- 
ysis of the 22 sets of acoustic focusing experiments and 
shown in the Supplemental Material. As in Ref. [l5| , the 
measured microparticle velocities u p are thus represented 
on a discrete x n x y m micro-PIV grid 



u p = u p (x n ,y m ) 



(19) 



All measured velocities presented in the following are nor- 
malized to their values at U pp = 1 V using the voltage- 
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FIG. 6. Experimental micro-PIV measurement of acousto- 
phoresis in the horizontal ultrasound symmetry plane of 1-um- 
diameter polystyrene particles suspended in water and driven 
at the piezo frequency / = 1.940 MHz corresponding to A/2 = 
w, temperature 25 °C, and voltage Up P = 7.94 V. The arrows 
represent the measured velocity vectors u p and the colors 
their magnitude it p normalized to Up P , see Eq. (|20|) . 



squared law [14[, 



1 V 

U*~ 
^ pp 



p* 



(20) 



where the asterisk denotes the actual measured values. 
As a result, the extracted velocity amplitudes and acous- 
tic energy densities are normalized as well, 



E„ 



1 V 

TJ* I ac- 
PP- 



(21a) 
(21b) 



The actual peak-to-peak values of the applied voltage 
C/*p for all four experimental series are given in the Sup- 
plemental Material. 



A. Excitation of a ID transverse standing wave 

In Figs.[5]and[7]we verify experimentally that the acou- 
stophoretic particle velocity is of the predicted sinusoidal 
form given in Eq. (|12p and resulting from a ID trans- 
verse standing wave. For the actual applied voltage of 
U* = 7.94 V the maximum velocity was measured to 



be 1.77 mms , which by Eq. (|20|) is normalized to the 
maximum velocity itP* ax = 28 urns -1 seen in Fig. [S] 

A detailed analysis of the measured velocity field re- 
veals three main points: (i) The average of the ratio of 
the axial to the transverse velocity component is practi- 
cally zero, (|mP/uP|) < 5 %, (ii) the maximum particle 
velocity along any line with a given axial grid point co- 
ordinate x m varies less than 6 % as a function of x m , 
and (iii) the axial average (u p ) x of the transverse veloc- 
ity component u p is well fitted within small errorbars 
(< 1 %) by Eq. CE2D. 



8 



20 ■ 


(a) 




1 , 






^ *w particles 













20- 




/ = 1.940 MHz ^^- w * arX ^ 





=1 



500 


(b) 








f particles 




o- 








500 




/ = 1.940 MHz ^^^T 





0.05 


(c) 


^ X particles 




o- 


7 


/ = 3.900%. / \ J 




0.05 ■ 




MHz 





-200 



-100 





j/ \fj,m] 



100 



200 



FIG. 7. Measured average {uy) x (black and red points) and 
its standard deviation a({uy) x ) (error bars) of the trans- 
verse velocity Uy of polystyrene particles in water. The si- 
nusoidal prediction Eq. (|12[) f° r w p (blue line) is fitted to 
data points {u^) x away from the side walls (red points), (a) 
Streaming-dominated acoustophoresis for 2a = 1 um driven 
at / = 1.940 MHz (A/2 = w, same data as in Fig. [6]). 
(b) Radiation-dominated acoustophoresis for 2a = 10 um 
at / = 1.940 MHz (A/2 = w). (c) Streaming-dominated 
acoustophoresis for 2a = 1 um at / = 3.900 MHz (A = w). 



B. Measuring the velocity amplitude 

In Fig.[7Ua) we plot the axial average (vR) x of the trans- 
verse velocity component (black and red points) and 
its standard deviation a((u^} x ) (error bars) for the ve- 
locity field shown in Fig. [S] at the standing half- wave res- 
onance frequency / = 1.940 MHz for the 1-um-diamcter 
streaming-dominated particles (series MQ1). The mea- 
sured velocities away from the side walls (red points) 
are fitted well by the predicted sinusoidal velocity pro- 
file u a sm(2ky) (blue curve) Eq. (fT2|) for fixed wavelength 
A = 2ir/k = 2w and using u a as the only fitting param- 
eter. Velocities close to the side walls (black points) are 
discarded due to their interaction with the side walls. As 
seen numerically in Fig. HI the no-slip boundary condition 
on the side walls of the rectangular geometry suppresses 
the streaming velocity near the side walls relative to si- 
nusoidal velocity profile of the parallel-plate geometry. 

As shown in Fig. EJb), the theoretical prediction also 
fits well the measured velocities away from the side walls 
for the large radiation-dominated 10-um-diameter parti- 
cles (series MQ1, A/2 = w). Likewise, as seen in Fig.jTJc), 
a good fit is also obtained for the 1-um-diamcter parti- 
cles away from the side walls at the standing full-wave 
frequency / = 3.900 MHz (series MQ2, A = w). 

Given this strong support for the presence of stand- 
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FIG. 8. Measured and normalized acoustophoretic velocity 
amplitudes u a /ito as function of particle diameter 2a (colored 
symbols) and wall-drag-corrected particle diameter 2a w d 
(gray symbols) for the four experiment series (a) MQ0, 
(b) MQl, (c) MQ2 and (d) Gil described in Section ITTTBl 
The characteristic velocity amplitude u is determined from 
fitting Eq. ()12p to the gray points in each series using Z? ac 
and s as fitting parameters resulting in the values listed in 
Table lIVT a - ! (green dashed lines). The standard deviation 
of u a /ti is less than the symbol size and the standard de- 
viations on the particle diameters are indicated as black lines. 



ing transverse waves, we use this standing-wave fitting 
procedure to determine the velocity amplitude u a in the 
following analysis of the acoustophoretic particle velocity. 

In spite of the normalization to the same driving volt- 
age of 1 V, the velocity amplitude of the half-wave res- 
onance in Fig. [TJa) is 400 times larger than that of the 
full- wave resonance in Fig. [TJc). This is due to a differ- 
ence in coupling to the piezo and in dissipation. 



C. Velocity as function of particle diameter 

To analyze in detail the transverse velocity amplitude 
it a in all four series MQ0, MQl, MQ2, and Gil, we return 
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TABLE IV. Measured acoustic energy densities E ac normal- 
ized to U pp = 1 V and streaming coefficient s. 



(a) Un-weighted fit to all points, see Fig 


El 


Susp., freq. 


S ac [J m- 3 ] 


s 


MQO, 1.940 MHz 


52.306 ± 0.918 


0.222 ± 0.025 


MQ1, 1.940 MHz 


31.807 ± 0.569 


0.247 ± 0.071 


MQ2, 3.900 MHz 


0.070 ± 0.001 


0.262 ± 0.125 


Gil, 2.027 MHz 


2.420 ± 0.020 


0.184 ± 0.012 


(b) Based on particles with 2a = 0.6 um 


and 2a = 10 jim 


Susp., freq. 


£ ac [J m- 3 £ 




MQ1, 1.940 MHz 


32.436 ± 1.282 


0.182 ± 0.008 


MQ2, 3.900 MHz 


0.071 ± 0.003 


0.205 ± 0.008 


Gil, 2.027 MHz 


2.559 ± 0.110 


0.186 ± 0.008 



a Eq. J23J 
b Eq. |QD 



to the wall-enhanced drag coefficient x of Section Hi Dl In 
general, x depends in a non-linear way on the motion and 
position of the particle relative to the rigid walls. How- 
ever, in Section IIII CI we established that the majority 
of the observed particles reside in the middle half of the 
channel, and in our standing-wave fitting procedure for 
u a in Section IIV Bl we discarded particles close to the 
side walls. Consequently, given this and the values of \ 
in Table U it is a good approximation to assume that all 
involved particles have the same wall correction factor, 
namely the symmetry-plane, parallel-motion factor, 



X 



paral 
Xz=0 • 



(22) 



As the drag-correction only enters on the radiation- 
induced term in Eq. (|15p. we introduce a wall-drag- 
corrected particle size a wc i = (Xz=o') _ ^ a - 

To determine the acoustic energy density E ac and the 
streaming coefficient s we plot in Fig. [8j for each of the 
four experiment series, u a /u versus the particle diame- 
ter 2a (colored symbols) and wall-drag-corrected particle 
diameter 2a w( j (gray symbols). The characteristic veloc- 
ity amplitude u is determined in each scries by fitting 
the wall-drag-corrected data points to Eq. (fT2"]) using E ac 
and s as fitting parameters. In all four experiment se- 
ries, a clear a 2 -dependence is seen. Notice further that 
the velocities follow almost the same distribution around 
the fitted line in all series. This we suspect may be due 
to systematic errors, e.g. that the 5-pm-diamcter parti- 
cles are slightly underestimated (see the Coulter data in 
Supplemental Material). The resulting fitting parameters 
E ac and s are listed in Tabic |lVja). The energy densities 
normalized to U pp = 1 V, see Eq. (|21b|) varies with more 
than a factor 700 due to a large difference in the strength 
of the excited resonances. According to the predictions in 
Section |TI] the streaming coefficient s should be constant, 
but experimentally it varies from 0.18 to 0.25. However, 
taking the fitting uncertainties into account in a weighted 
average, leads to (s) w = 0.192±0.010 close to sf « 0.194 




O Gil, 2.027 MHz 
□ MQ1, 1.940 MHz 
V MQ2, 3.900 MHz 

(a wd /a ) 2 + sf 
std. dev. 



10 15 
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FIG. 9. Normalized acoustophoretic particle velocities u a /u 
versus normalized particle size squared (o w d/oo) 2 . 



of Eq. P]). 

Another approach for extracting E ac and s is to as- 
sume that the smallest particles 2a = 0.6 jim are in- 
fluenced only by the streaming-induced drag. If so, the 
velocity of the largest 2a = 10 jim particle has a stream- 
ing component of less than 6 %, see the measured ratios 
u a 6 lim / M i° pm m Table [V] Therefore, we further assume 
that the 10 um-diamcter particles are influenced solely 
by the radiation force, and from u T a d = u a (a/a ) 2 we 
determine the acoustic energy density as 



3?7c 



.10 pm 



2 $w 



(23) 



Knowing the acoustic energy density, we use Eq. ([7]) to 
calculate the streaming coefficient s from u s a r = u s as 



_ Po c o 0.6 um 
4£ " a ■ 



(24) 



Assuming that the largest error is due to the dispersion in 
particle size, we obtain the results listed in Table ITvT b). 
The acoustic energy densities are close to the ones ex- 
tracted from the fits in Fig. [S] and the geometric stream- 
ing coefficient varies from 0.180 to 0.203 with an weighted 
average of (s) w = 0.191±0.005. Note that using Eqs. ([23]) 
and ((2"4"]) , we only need to consider the dispersion of the 
10-pm-diameter particles, which results in a more reliable 
estimate of s. 

We use the acoustic energy densities in Tabic HVf b) 
together with the material parameters in Table |TT] to cal- 
culate u and ao, Eq. ([SJ, for each of the experiment 
series MQ1, MQ2, and Gil. According to the theoreti- 
cal prediction in Eq. (|12p . all data points must fall on a 
straight line of unity slope and intersection s if plotted 
as the normalized velocity amplitude u a /u as function 
of the normalized particle radius squared (a/ao) 2 . The 



10 



TABLE V. Relative particle velocities. 



Susp., freq. 


„0.6 B m 


„10 w 




10 pm 


ff.e Jim 


1U um (J.fi um 


MQ1, 1.940 MHz 


0.020 


49.6 


0.294 


MQ2, 3.900 MHz 


0.011 


88.4 


0.291 


Gil, 2.027 MHz 


0.061 


16.4 


0.270 



plot is shown in Fig. [5] showing good agreement with the 
theoretical prediction using « 0.194. 

D. Velocity ratios 

In Table [V] we list velocity ratios for different particle 
sizes in the experiment series MQ1, MQ2, and Gil. 

From Eq. (|12|) we expect u a — w| tr cx a 2 leading to 
the prediction (u 5 a vm - uf T ) / (n™ * m -uf T ) = (5/10) 2 = 
0.25. If we assume that the smallest 0.6-um-diameter 
particles are only influenced by the acoustic stream- 
ing, we have u| tr = 6 11111 . We can therefore test 
the just mentioned hypothesis by calculating (u^ 11111 — 

M 0.6 »m)/( u 10 pm _ ^0.6 V my The results are ligted m th(J 

third column in Table |V1 where we obtain values ranging 
from 0.27 to to 0.29, or a deviation of 8 to 18 %. 

Assuming that the smallest 0.6-pm-diameter parti- 
cles and the largest 10-um-diameter particles arc influ- 
enced only by the acoustic streaming and the acous- 
tic radiation force, respectively, we can estimate the 
ratio of radiation- and streaming-induced velocities as 
< ad /< tr = ul° lim /u°' 6 ym , which are listed in the sec- 
ond column in Table [V] First, we notice that the ratio 
increases by a factor of 88.4/49.6 = 1.8 as we increase the 
frequency by a factor of 3.900/1.940 = 2.0. This agrees 
well with a linear increase with frequency as predicted by 
Eq. (fTS"]) . Secondly, we notice that the ratio increases by 
a factor of 48.6/16.4 = 3.0 as we change the suspending 
medium from a 0.75:0.25 mixture of water and glycerol 
to pure water. According to Eq. ([T3]) M lad /w| tr increases 
linearly with ^/$ and from Table [IT] we obtain a pre- 
dicted ratio increase of 16.8/5.25 = 3.2, which matches 
well with the experimentally-estimated ratio. 

With these results we have gained experimental sup- 
port for the theoretical prediction of the velocity ratio 
given in Eq. (|13[) . 

V. DISCUSSION 

Our results verify experimentally the theoreti- 
cally predicted dependence of the magnitude of the 
acoustophoretic velocity in a microchannel on the viscos- 
ity of the suspending liquid, the acoustic contrast factor 

and the ultrasound frequency. 

For most situations involving cells, isotonic solutions 
are used such as PBS, sodium chloride, or blood plasma. 
Direct manipulation of small particles such as bacteria in 



plasma or in buffers of high levels of protein is problem- 
atic, primarily due to the high viscosity of those media. 
When possible, these media should be exchanged prior 
to manipulation to increase the potential for success. 

Increasing the frequency in the system would allow for 
a reduction of the critical diameter of particles. One 
adverse effect of a higher frequency is that the channel 
width must be narrowed down which affect the through- 
put in the system. This problem can, however, be over- 
come by designing a channel of high aspect ratio, where 
the resonance occurs over the smallest dimension [39j . 

Another benefit of high aspect ratio channels was 
pointed out by Muller et al. (2J]. Since the acoustic 
streaming emanates from the walls perpendicular to the 
wave propagation, here the top and bottom, a high chan- 
nel leads to a weaker average streaming field in the center. 

The measurements of particle velocities for polystyrene 
particles ranging from 0.6 um to 10 pm give no support 
to previous measurements presented by Yasuda and Ka- 
makura in 1997. Their rather spectacular result was 
that particles below a certain size move faster than larger 
particles do. From the experiments reported herein it is 
clear, however, that the motion of particles indeed can be 
well described with the analyses presented by Rayleigh 
[23l | , Yosioka and Kawasima [TtJ , and Gorkov [18| . 

The uncertainties in the measured particle velocities 
may in particular be due to the following four causes: 

(i) Variations in particle density p p and compressibil- 
ity K p as function of particle producer (Fluka and G. 
Kisker) or batch, (ii) deviations from normal distributed 
particle sizes as shown in Supplemental Material Fig. 1, 

(ii) local fluctuations in the bead concentrations leading 
to particle-particle interactions, and (iv) viscosity varia- 
tions induced by the suspended particles. 

To better understand the nature of acoustic streaming 
in microchannels the streaming field should be mapped 
for the channel cross section and along the whole length 
of the channel. As reported by Hagsater et al. [4l[ and 
Augustsson et al. [l5[ the acoustic field can deviate dra- 
matically from the very uniform one dimensional field 
reported herein. Non symmetrical acoustic fields can be 
expected to generate far more complex streaming fields. 



VI. CONCLUSIONS 

We have investigated the motion of microparticles due 
to acoustic radiation and acoustic streaming inside a 
liquid-filled long, straight rectangular channel of width 
w and height h driven by an ultrasound standing wave 
of wavelength A . 

Fortuitously, the simple analytical expression derived 
by Lord Rayleigh for the streaming velocity in an isother- 
mal liquid slab between two infinite parallel plates ful- 
filling A > ft, is a good approximation for the specific 
rectangular channel of Fig. [1] containing a thermoviscous 
liquid and fulfilling A ~ h. The reduction in velocity 
obtained when substituting the parallel plates with the 
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rectangular geometry is almost perfectly compensated 
for by the enhancement in velocity from substituting the 
isothermal liquid by the thermoviscous one. 

A theoretical prediction was made, Eq. (|12p , for the de- 
pendence of the radiation- and streaming-induced veloci- 
ties on the size of the particles, the ultrasound frequency, 
the viscosity of the suspending liquid, and the acoustic 
contrast factor. This prediction was found to be in ex- 
cellent agreement with experimental findings as shown 
by the collapse after re-scaling of data from 22 differ- 
ent measurement on the same line in Fig. [9j The results 
have bearing on acoustophoretic manipulation strategies 
for sub-micrometer biological particles such as bacteria 
and vira, which are too small to be handled using the 
present manifestation of this technique. We can con- 
clude that increasing the ultrasound frequency, increase 



of the channel aspect ratio, and lowering the viscosity of 
the suspending fluid is probably the most viable route to 
conduct such manipulation. 
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